Nonlocal fluctuations and control of dimer entanglement dynamics 



Cristian E. Susa and John H. ReinaQ 

Departamento de Ftsica, Universidad del Valle, A. A. 25360, Call, Colombia 
(Dated: February 8, 2010) 

We report on the dissipative dynamics of an entangled, bipartite interacting system. We show 
how to induce and control the so-called early stage disentanglement (and the 'delayed' entanglement 
generation) dynamics by means of a driving laser field. We demonstrate that some of the features 
currently associated with pure non-Markovian effects in such entanglement behavior can actually 
take place in Markovian environments if background noise QED fluctuations are considered. We 
illustrate this for the case of a dimer interacting molecular system for which emission rates, inter- 
action strength, and radiative corrections have been previously measured. We also show that even 
in the absence of collective decay mechanisms and qubit-qubit interactions, the entanglement still 
exhibits collapse-revival behavior. Our results indicate that zero point energy fluctuations should 
be taken into account when formulating precise entanglement dynamics statements. 
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I. INTRODUCTION 

Entanglement is arguably the most striking feature of 
quantum phenomena [1 j, and is a crucial physical re- 
source to quantum computing [2] , and quantum commu- 
nication via protocols such as teleportation [3] and cryp- 
tography [4 . Works on entanglement dynamics have re- 
cently shown that some initial bipartite entangled states 
can decay to zero in a finite time much shorter than that 
of their spontaneous emission jSUTTj . This phenomenon- 
early stage disentanglement (ESD)- is a quantum feature 
that signals an unusual dissipative dynamics of purely 
non-local quantum correlations [TJ [3] . It has also re- 
cently been shown that, in specific scenarios, it is possible 
to produce a delayed creation of entanglement, the so- 
called entanglement sudden birth (ESB) [H[T3]. These 
entanglement features have been explored in both the 
Markovian [14] and Non-Markovian [15] regimes in sev- 
eral different contexts [TTHT6] . Even more recently, there 
have been reports on the quantification of non-classical 
correlations of bipartite qubit systems in interaction with 
non-Markovian reservoirs [17 by means of the quantum 
discord [18] dynamics, in contrast to that of entanglement 
evolution [T7] . 

In this work, we show how the ESD and ESB phe- 
nomena can be externally coherently controlled with a 
laser field. We also show that the entanglement dynamics 
is permanently affected by the background noise due to 
nonlocal vacuum fluctuations, a fact that permanently af- 
fects the ESD dynamical profile. Our results give a prac- 
tical prescription, within the Born-Markov formalism, for 
the identification and control of ESD, ESB, and collapse- 
revival phenomena. This shows that previously identified 
entanglement dynamics associated to non-Markovian be- 
havior [15], can actually take place within a Markovian 
dynamics description. 
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Most of the interacting two-qubit systems can be de- 
scribed by a Hamiltonian of the type [19J Hs — ^0+^12, 
where Ho denotes the free particle term, and H12 = 
h(J x a x ^ (g) ai 2) + JyGy^ ® cr^ + J Z (J^ ® cr^ 2) ) , writ- 
ten in terms of the Pauli matrices, denotes a generic, 
anisotropic qubit-qubit interaction, and h is Planck's 
constant. In many physical processes, such as energy 
transfer in biomolecular systems [16] [20] or interacting 
artificial atoms [21] . to name but a few, J x = J y = J xyi 
and the strength ratio J = J xy / J z sets the type of inter- 
action experienced by the coupled qubits, and defines the 
type of entanglement that can be 'naturally' [19] gener- 
ated in the bipartite system. 

In this article, we consider the case of dipole-dipole 
(d-d) interacting single molecules of transition frequen- 
cies V{ (molecule i) in interaction with the quantized ra- 
diation field, and externally driven by a coherent laser 
field [221 [23]. We denote by |0 f >, and |1;), i = 1,2, 
the ground and excited state of molecule z, respectively, 

and hence Ho = — \v\(j)t^ — \v2(?^ • The two two- level 
molecules are separated by the vector 1*12 and are char- 
acterized by transition dipole moments jli = (0^ 1 1^) , 
with dipole operators D^, and spontaneous emission rates 
Ti . The system interaction Hamiltonian can be written in 
the computational basis of direct product states \i) <g> \j) 
= 0,1) as Hs = Hq + -Hi 2? where H12 is set by 
the dipole coupled molecules of interaction energy hV\2 

(J x = Jy = V 12 , J Z =0). 



II. DIMER DISSIPATIVE DYNAMICS 

Evidence of the molecules dipolar coupling, and the 
generation of sub- and super-radiant states, as well as 
an effective shift of the doubly excited state, have been 
experimentally reported in [22]. The latter has been ac- 
counted for in H\2 as a phenomenological, effective A e - 
shift of the doubly-excited state |11) = 1 1 1 ) (g) 1 12) ? due 
to the vacuum fluctuations [23]. The dimer is externally 
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controlled by a coherent driving field which acts on each 
of the molecules with a coupling amplitude Mi = —/j,i-Ei, 
and a frequency ujl- The light-matter interaction Hamil- 



tonian# L = ^(^e 



(i) Aoj L t _ 



-afe 



where \±i is the 



i-th transition dipole moment and Ei is the amplitude of 
the coherent driving acting on molecule i (at Vi). 

The dimer dissipative dynamics is appropriately de- 
scribed within the Born-Markov formalism, with H = 
Hs + Hl, by the quantum master equation 



[H,p]+L(p), 



(1) 



with dissipative Lindblad super-operator given by 



m = -^K'^+^^p-^V?') (2) 
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1^21 (2) (1) . (2) (1) A Q (2) A 
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where 1^ = I\ = nA 12 -i/2 = ncjf ||^|| 2 /(3e ^c 3 ) are the 
spontaneous emission rates (A\2\U i = 1,2 are the Ein- 
stein vacuum coefficients), n is the refraction index of the 
dispersive medium in which the molecules are embed- 
ded (for example, a paraterphenyl crystal in [22J), and 
Ti2 = 1^21 are the collective spontaneous emission rates, 
arising from the coupling between the molecules through 
the vacuum field [24 ; af = |1»><0*|, and a® = |0*><1*| 
are the raising and lowering Pauli operators acting on 
molecule i. 

In particular, in the near field approximation (7*12 <C 
Xl) the d-d interaction energy is given by hV\2 = 

3h £W~ 2 [Ai • A2 - 3(Ai • fi 2 )(A2 • fi 2 )], and T 12 = 
yTiT^Ai ' A2 : the maximum (minimum) interaction 
strengths V12 are obtained for parallel (perpendicular) 
dipole moments [33]. z = nk^r^] fii, and 1*12 are unit 
vectors along the transition dipole moments and along 
the relative separation between the molecules. 

We stress that Eq. ([I]) is more general than that used 
in [TTJ [15] in the sense that it accounts for qubit-qubit 
coupling, coherent external driving, and nonlocal fluctu- 
ations due to radiative corrections 24 . 



III. RESONANCE FLUORESCENCE 

To evaluate the system's physical observables, we per- 
form the calculation of the occupation probabilities in 
the steady state. We compute this for several different 
experimental conditions realizable in the laboratory, as 
shown in Fig. 1. These involve weak and strong coherent 
laser driving (e.g., Fig. 1(c)), different dipole interaction 
strengths (Figs. 1(a), and 1(d)), and different dimer 's 



resonance frequencies (Fig. 1(b)). In particular, we ac- 
count for the fluorescence spectrum measured for the pair 
of coupled molecules reported in Ref. [22 . 

Our calculations are performed in the standard two- 
qubit computational basis {|00), |01), |10), |11)} = 
\h) ® b'2))? an d the corresponding dimer density matrix 
elements are denoted as Pij,kU with associated popula- 
tions pij^j (density matrix elements i,j = 0,1. 
The calculations here presented are amenable to experi- 
mental verification in a dimer molecular set-up where the 
task of optically resolving the single molecules by means 
of laser spectroscopy techniques is possible. This kind of 
experiment has been performed, for instance, by Hettich 
et al. [22], and the detection and measurement of cou- 
pled single terrylene molecules embedded in an organic 
(para-terphenyl) crystal (thickness of about 250 nm) has 
been reported. There the detection was carried out by 
means of fluorescence excitation spectroscopy at a tem- 
perature ~ 1.4 K. Such a recorded fluorescence spectra, 
at high excitation laser intensity, exhibits the appearance 
of an 'additional' peak between those associated to the 
individual transition frequencies and V2, which is the 
experimental evidence of the dipole-dipole coupling V12 
between the qubits. 

The continuous curves shown in Fig. [TJa) are in perfect 
agreement with the experimental results reported in [22] , 
for the same set of implemented physical parameters: 
difference in dipole transition frequencies A_ = 2320 
MHz, dipolar coupling strength V12 = 950 MHz, indi- 
vidual and collective decay rates Ti = T 2 = 27r x 50 
MHz, Ti2 = 2tt x 9 MHz; with A e = -160 MHz, for a 
laser coupling ^ = 200 MHz. Here, we have calculated, 
directly from Eq. ([T]), the steady state occupation prob- 
abilities for the states |01), |10), and |11) as a function of 
the laser detuning. We can also see from Fig. [TJa), the 
effect due to a different dipolar coupling, as shown by the 
dashed curves, for V12 = 50 MHz. Since this implies a 
change in the relative position of the coupled dipoles, this 
is reflected in a shift of the resonant peaks at frequencies 
^1, and V2, and a decreasing of the signal intensity com- 
ing from the two-photon emission process (intermediate 
peak) due to the weaker dipolar coupling. 

The case of molecules with equal transition frequen- 
cies (A_ = 0) is illustrated by the continuous curves in 
Fig. [jjb), for V12 = 950 MHz. As expected, only two 
peaks appear: one associated to the resonance condition 
v\ = z/2, and the intermediate one, associated to the dipo- 
lar coupling. By bringing the dipoles into resonance, a 
clear increment in the two-photon resonance signal, in 
comparison to the experimentally measured spectra [22] 
(calculated as the dashed curve in Fig. [ljb)) is observed, 
thus indicating a stronger dipole-dipole coupling. 

The resonance fluorescence spectrum reported in [22] 
can be theoretically reproduced by plotting the steady 
state quantity poi,oi + Pio,io + 2pn j n, as shown in Fig. 
[TJc). This shows the appearance of the two-photon reso- 
nance as the laser excitation intensity is increased. This 
peak is shifted from the laser detuning A + = due to the 
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FIG. 1: Steady state occupation probabilities for the dimer system. The detunings A_ = v\ — z/ 2 , and A+/2 = (y\ + zy 2 )/2 — v^. 
a) Populations pij,ij for two different dipolar couplings, A_ = 2320 MHz. b) Probabilities for a fixed dipolar coupling, = 950 
MHz, and different molecular resonance conditions. In both cases, ri = Y2 = 2tt x 50 MHz, Y12 = 2tt x 9 MHz, A e = —160 MHz 
and ti — 200 MHz. c) Fluorescence spectrum given in terms of the occupation probabilities, as a function of the laser coupling 
strength and the laser detuning A+/2. In MHz, Y% = 187T, Y12 = 97T, A_ = 2320, and A e = —160. d) Zoom of the two photon 
(intermediate) resonance for £i — 100 MHz. 



the vacuum fluctuation noise, an effective shift that has 
been experimentally [22] and theoretically [23] estimated 
at around —160 MHz. The calculations here presented 
have been done for the independent and collective de- 
cay rates r, = 18tt MHz, r i2 = 9tt MHz. Figure [ltd) 
shows a zoom of the intermediate resonance of Fig. (TTc) , 
for a laser coupling ti = 100 MHz, and different dipolar 
couplings. This shows a clear increment in the size of 
the resonance as the qubit coupling is increased. If this 
graph were to be calculated for ii = 10 MHz (not shown), 
the probability of simultaneous excitation pn ? n appears 
multiplied by the scale factor 10 -5 . Hence, the resonance 
signal can be tailored according to designed structural 
molecular parameters. The difference in the half-width 
of the reported peaks has its origin in the constructive 
and destructive interference of the different possible re- 
laxation pathways [25] . 



IV. DIMER ENTANGLEMENT EVOLUTION 

As a norm for the entanglement quantification of the 
dipole-coupled molecules, we use Wootter's concurrence 
C(p) = max{0, Ai — A2 — A3 — A4}, where the A^'s are 
the square roots of the eigenvalues of the non-hermitian 



matrix pp, arranged in decreasing order, with p = a y 



(!) , 



matrix 



(2) 

(j y ' [26] . We consider the initial state density 



p(0) = 



a w 

b z 

z* c 

\w* d 



(3) 



as this allows the initialization of a broad class of en- 
tangled states [HJE]. In the standard two-qubit com- 
putational basis, the maximally entangled Bell state 
|^+) ^(|01) + 1 10>) is obtained for a = w = d = 0; 
b = c = z = 1/2. All the other states of the Bell basis, 
as well as the Werner mixed state [14] can be produced 
in a similar manner. 



A. Coherent control dynamics 

We consider the dynamics of the initial entangled state 
\i/jo(a)) = y/a\01) + v^|10)? where a is a real number, 
VP — \\fP\ e ^ is a complex number of phase 0, and 
a + |/3 1 = 1. This choice follows the identification of 
the eigenstates of the system's Hamiltonian Hs, in the 
absence of external driving. If no cooperative phenomena 
due to V12 (decay rate Yi 2 ) take place (for example, due 
to a long separation between the molecules, V12 ~ 0), 
in the absence of laser driving, the entanglement of this 
type of state has a natural asymptotic time decay due to 
the emission rates Ya = Yi of the individual molecules. 

The numerical solution of the master equation ([!]) al- 
lows the quantification of the dimer entanglement dy- 
namics, as shown in Fig. [2] where the external field 
has been used to produce a coherent control of entangle- 
ment. Figure [2^a) shows that for a strong laser coupling 
{ii — 500 MHz) the concurrence abruptly decays to zero 
for a time less than tp, where tr = 1/Y is the time asso- 
ciated to the decay rates Yi = Y for all the initial state 
configurations (all a). In this work, A + , A_, A e , V12, 
and £i are given in frequency (y = uj/2tt) units, whereas 
and Yij appear divided by a 27r-factor to consistently give 




FIG. 2: Laser-driven entanglement dynamics control, a) 
Early stage disentanglement: A+ = 0, and £\ — £ 2 = 
500 MHz. b) The molecules couple to the laser with different 
strengths: h = 300 MHz, and £ 2 = 500 MHz. Figures c) 
A+ = A_, and d) A+ = — A_, show the difference in entan- 
glement behavior when the laser resonant condition is varied, 
for £ x = 300 MHz, £ 2 = 500 MHz. A e = 0. In all the graphs, 
r x = T 2 = 2tt x 50 MHz, Ti2 = 2tt x 9 MHz, A_ = 2320 MHz, 
and V12 = 950 MHz. 



currence dynamics-the sudden birth of entanglement. In 
these graphs, the system is prepared in the initial sepa- 
rable state 1 1 1 ) g) (y/a\02) + \//?|l2)), ano ^ after a certain 
time, which depends on the type of input superposition 
10 (i.e., on the value of a), entanglement appears in the 
system (see Figs. [3^a), and (c)). Note that, due to the 
detuning between the molecular frequencies, the different 
resonant laser tuning of the molecules produces different 
entanglement behavior. 

Figure [3^a) shows a delayed ESB for a = and de- 
values close to zero. The fast oscillations of the concur- 
rence in Fig. [3ja) are plotted in its inset for about 1/16 of 
the total time-frame considered in the main figure. The 
graphic) shows the effect due to the change in the laser 
i.oresonance molecular tuning. Here, A + = — A_, and an 
appreciable 'a-region' of early stage disentanglement is 
suppressed, accompanied by an increment in the amount 
of generated stationary entanglement. As a major fea- 
ture, we notice that due to the laser addressing strategy 
(let's recall that the two dipoles have different transition 
energies), the initial states |0) <8> and |1) g) |0) (a = 1) 
are more prone to a quick generation of entangled states 
than the initial condition |1) (g) |1) (a = 0); in the latter 
there is an initial lack of entanglement generation, and 
ESB behavior takes place. In the former, the quick en- 
tanglement generation can be understood from the fact 
that a coherent superposition of the states |0) g) |1) and 
|1) g) |0) is naturally supported by the dipole-dipole in- 
teraction; in fact, in the absence of laser tuning, such a 
superposition of states produce some of the eigenstates of 
the system's Hamiltonian. The time scale for which ESB 



the same units. This early stage disentanglement is due 
to the vacuum fluctuations affecting each molecule sepa- 
rately, and radically changes when the coupling between 
the laser and the molecules £\ 7^ £2, as shown in Fig. 
f2^b) . There, for all a, a finite, stationary entanglement 
C(p) remains for long times. This behavior occurs af-C(p) ; 
ter a period for which the system experiences a complete 
disentanglement — C(p) = for all a, and can be further 
controlled, as shown in Figs. [2|c) and (d), by tuning the 
molecule resonance conditions with the laser field. 

In Fig. [2|c), the ESD behavior is partially suppressed 
for finite a (0.8 < a < 1), when the laser is on reso- 
nance with the qubit 2 (vl — v<i), that is, for A + = A_ 



(y\ > 1/2). In this case, the asymptotic limit of the con- 
currence is increased to ~ 0.25. In Fig |2|d), the 
quency of the external field has been tuned to be in res- 
onance with that of qubit 1 (A + = — A_): the ESD 
behavior is now completely suppressed for all a, and the 
entanglement shows an oscillatory decay. For long times, 
and for all a, the concurrence approaches the same sta- 
tionary entanglement as before. The fast oscillations of 
the concurrence are due to the competition between both 
the dipolar coupling and the laser-molecule interactions. 
Thus, the results of Figs. [2|c) and (d) give a mechanism 
for coherent control of the dimer's ESD behavior, and for 
the generation of steady-state entanglement. 

Figure [3] shows another interesting behavior in the con- 
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FIG. 3: Entanglement dynamics and ESB. Initial separable 
state |1}(a/^|0) + v^|l», (a) A+ = A_, h = £ 2 = 100 MHz; 
(b) A+ = A_, £ x = 300 MHz, £ 2 500 MHz. (c) A+ 
-A_, £i = 100 MHz; (d) A+ = -A_, h = 300 MHz, and 
£ 2 = 500 MHz. ri = T 2 = 100tt MHz, T 12 = 18tt MHz, 
A_ = 2320 MHz, A e = -160 MHz, and V12 = 950 MHz. 
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C(p) 




FIG. 4: Entanglement dynamics for initial entangled states 
given by Eq. (J3|, with a H> |, b = c = z = |, w = 0, 
and d = under spontaneous emission. A_ = 2320 MHz, 
A+/2 = (yx + zy 2 )/2 « 1 x 10 9 MHz (optical transition). The 
decay rates Y% = V = 1007T MHz, T12 = 0; the remaining 
parameters (V12, and A e ) are set to zero. 



takes place depends on the qubit resonance conditions: 
for the initial condition |1) <g> |1), the ESB appearance in 
Fig. |3^a) takes a longer time than that calculated for 
Fig. |gc). 

Figures |3jb) and (d) show the effect due to increas- 
ing the laser intensity. A quicker generation of entan- 
glement at around a = is now produced (cf. Figs. 
3 
3 



a) and (c)), with a subsequent generation of ESD (Fig. 

b) ), and a partial ESD suppression (Fig. |3jd)), fol- 
lowed by an overall increment in the amount of gener- 
ated stationary entanglement, C ~ 0.2. These graphs 
manifest the effect due to the tuning of the molecu- 
lar resonance: if A + = — A_ (Fig. |3jd)) the time- 
frame for the region of zero concurrence of Fig. [3jc) 
is almost completely suppressed, in addition of a larger 
number of entanglement oscillations before the steady- 
state is reached. For the chosen initial condition, the 
A + = — A_-resonance condition favors the entanglement 
production and inhibits ESD behavior. The same con- 
clusion can be reached if we instead consider, as initial 
condition, the state (y/a\0) + v7?|l)) (g) |1) (not shown). 
The changes in the entanglement features shown in Fig. [3] 
are not only due to the chosen initial conditions and reso- 
nance features, but also to the fact that, in the considered 
scenarios, we have taken into account the effective energy 
shift A e = —160 MHz, which is associated to the dou- 
bly excited state, as reported in [22, 23]. The numerical 
values for the strength of the molecular interaction, the 
decay rates, and the molecular transition energies used 
in the simulations of Fig. [2] and [3] are those reported in 
the experiment of Ref. [22] . 



B. Dependence on the initial state preparation 

Further to our previous analysis, we plot in Fig. [4] the 
entanglement dynamics for an initial state of the form 



Eq. with a 4 |, b = c = z = 1/3; w = 0, and 

d = Due to the structure of the density matrix, 

the concurrence allows an analytical calculation of the 
observed ESD behavior, and the associated ESD time t 
(for each a) can be exactly estimated. For example, for 



0, t 



jUn 



2+V2 
2 



which means t ~ 1.7 ns < T 



our molecular set-up. This behavior is observed over a 
continuous range of a values, and for any 0<a<2/3is 
given by 



t 



V— 

r \2 



3a 



a + \fa? 



(4) 



For a > 2/3, ESD is suppressed and entanglement decays 
in an asymptotic fashion. Hence, the dissipative molecu- 
lar dynamics given by Eq. ([!]) exhibit a rich behavior and 
both an ESD dynamics and an asymptotic decay of en- 
tanglement can be induced in the dimer system by tuning 
the initial condition (a value). 

In Fig. [5] we consider the entanglement dynamics for 
the general initial product state of qubit superpositions 

1*0(7)) = (V7IO1) + Vs\u)) ® (x/CIOs) + Vv\h)), (5) 

where, as before, 7 + \S\ = 1, ( + |^| = 1. Here, we 
plot the corresponding concurrence as a function of 7, 
for C = 0, 1/2, 1 (vT > 0). We note that different initial 
preparations for the disentangled type of states lead to 
very different entanglement behavior. For instance, Fig. 
[5ja) exhibits a sudden birth of entanglement dynamics 
for low values of 7, whereas for their intermediate values 
there are collapses and revivals that eventually become 
suppressed as 7 is increased; from this point, coherent 
oscillations of entanglement persist to finally decay and 
then converge to a stationary value. In fact, the concur- 
rence evolves, for all 7, towards a finite non-zero station- 
ary value due to the inter-qubit coupling strength and the 
driving field. Also note that for short times the molec- 
ular entanglement reaches its maximum for 7 ^ 1. Fig. 
|5|b) illustrates the case of equal superpositions in both 
qubits. Unlike the previous case, here there is a finite 
entanglement oscillation at short times for 7 \-> 0, which 
can be understood in terms of the symmetry of the initial 
state. Here there is no presence of ESB behavior. Due to 
the initial condition, Fig. [5^c) shows, in contrast to Fig. 
|5|a), the highest entanglement oscillations for 7^0. In 
all cases, the oscillations approach the same stationary 
value, as previously observed. Notice that in Fig. [5] we 
have considered a non-zero value for V12 and for the A e - 
shift, in contrast to the graph of Fig. [4] The parameters 
used in this simulation have been taken from the exper- 
imental data of Ref. [22]. We have recently learned, 
although in a different context, that steady state entan- 
glement is also produced if counter-rotating wave terms 
are kept in the original system Hamiltonian [8j [27] . 

Figure [6] follows a similar analysis as before, but for 
the case in which there is no physical interaction between 
the molecules, and therefore no collective decay. We con- 
sider instead the effects due to the vacuum fluctuations 
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FIG. 5: Initial state preparation of product states ^0(7)) = (a/t|0> + VS\1)) ® (VC\0) + y/rjW), for a ) C = 0> b ) C = 1/2, c) 
C = 1. A_ = 2320 MHz, A+ = 0. The decay rates Yi = IOOtt MHz, Y 12 = 18tt MHz. h = 100 MHz, A e = -160 MHz, and 
V12 = 950 MHz. 
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FIG. 6: Initial product states ^0(7)) (Eq. |5])), in the absence of qubit-qubit interaction and collective decay (V12 = 0, Y12 = 0), 
for a) C = 0, b) C = 1/2, c) C = 1. A_ = 2320 MHz, A+ = 2638 MHz. Yi = 10tt MHz, h = 200 MHz, A e = -160 MHz. 



A e -the Lamb shift experienced by the molecular dou- 
bly excited state, and to the external driving field. The 
results are clearly different from those of Fig. [5j even 
though the qubits start from a separable state, they en- 
tangle due to the radiative corrections and the external 
field. The concurrence exhibits a similar behavior for the 
conditions a) £ = 0, and c) £ = 1, for which there is a 
maximum amount of generated entanglement, and ESD 
and collapse-revival behavior are present as a function 
of the initial qubit preparation. The scenario £ = 1/2 
(graph |6]^b)) does exhibit a different type of decay, and 
unlike the previous cases, there are not two "islands" of c ^ 
entanglement, and the oscillations persist the most for . 8 P 
71—^1. Notice that the time scale for such oscillations is ci.4w 
much longer than those of the graphs previously shown. 2 o! 



1. Laser coupling 

We address the entanglement effects due to the molec- 
ular coupling to the laser field, for the initial condition 
given by Eq. Figure [7] shows the entanglement evo- 
lution as a function of the laser field coupling strength £. 
In the figure, 7 = £ = 1/2, and the other parameters are 
given as in Fig [5} for a) V 12 = 950 MHz, and b) V 12 = 50 
MHz. Figure ma) shows that for the considered initial 
product state of superpositions there is an optimal laser 



coupling strength for which maximum stationary entan- 
glement is generated. Away from this optimal value, the 
stationary entanglement decreases monotonically until it 
reaches a certain intensity from which any entanglement 
correlation is washed away and the concurrence goes to 
zero after a few oscillations. This indicates that the laser 
can be used to tune and maximize stationary entangle- 
ment production, given an initial condition. 




FIG. 7: Concurrence dynamics as a function of the laser inten- 
sity ti = £, for the initial condition ^0(7)), with 7 = £ = 1/2. 
a) V12 950 MHz, b) V12 50 MHz. A_ 2320 MHz, 
A+ = 0, T = IOOtt MHz, Y 12 = 18tt MHz, A e = -160 MHz. 

In Fig. [7jb) we have decreased V12 to 50 MHz. In con- 
trast to the previous case, for low laser coupling strength, 
there are no entanglement oscillations before the con- 
currence reaches its minimum value: the concurrence 
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exhibits an island whose height depends on the laser 
strength. Unlike the previous case, here the maximum 
entanglement occurs at low intensity, due to the fact the 
qubit-qubit interaction is smaller than A e , which domi- 
nates the dynamics. The entanglement rapidly vanishes 
as the laser intensity increases, with the two-particle state 
becoming separable during a period whose length in- 
creases with the laser intensity. It then starts to build up 
towards a stationary degree of entanglement, for which 
there is also a maximum value which is directly depen- 
dent on the intensity of the laser. After this, the entan- 
glement gradually decays as the laser intensity is further 
increased (not shown) and eventually vanishes following 
a similar pattern to that of Fig. (7^a). 




FIG. 8: Entanglement dynamics as a function of the ini- 
tial condition and laser intensity, for the initial configuration 
\ipo(a)}. a) i = 200 MHz, b) a = 1/2. A_ = 2320 MHz, 
A+ = 2638 MHz, F = IOtt MHz, A e = -160 MHz. Viz = 
T12 = 0. 

For the sake of completeness, we next consider the sit- 
uation of non- interacting qubits (V12 = 0), and no collec- 
tive decay for the initial state \ipo(a)) = ^/a\01) + y//3\10) . 
Figure |8ja) shows that ESD and collapse-revival behav- 
ior take place in the entanglement evolution as the initial 
condition is varied, by means of tuning a (e.g., maximal 
entanglement is attained for a = 1/2). In contrast, Fig. 
[8^b) shows the dependence of entanglement on the laser 
intensity for the same initial state a = 1/2. Here, we 
distinguish three main dynamical scenarios. For low in- 
tensity values, there is an asymptotic decay of concur- 
rence. As the intensity increases, collapse-revivals occur 
whereas for larger values of field coupling, pure ESD phe- 
nomena take over the dynamics. The time scale plotted 
in Fig. M is to be compared with the 'natural' decay time 
tr = T . It is clearly shown that coherent laser driving 
can be used as a key ingredient for determining the major 
features of entanglement dynamics. 



C. Nonlocal fluctuations: Lamb shift correction 

We illustrate in Fig. [9] a major point of this work: 
in the scenario where the molecules are not subjected 
to the collective damping due to the vacuum field, dark 
periods and revival of entanglement appear. This result 



is in contrast to the revival of entanglement reported in 
[30] , where such a behavior is critically dependent on the 
collective decay. As is shown in our case, the revival phe- 
nomenon is strictly induced by the external laser control 
and the quantum noise fluctuations given by the Lamb 
shift 23]. 
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FIG. 9: Concurrence evolution as a function of the initial 
condition \^o(a)}. The qubits are far apart from each other 
and do not experience a collective decay (V12 = 0, T\2 = 0). 
They, however, due to the vacuum noise, experience A e = 
-160 MHz [22]. Ti = 4tt MHz, A_ = 2320 MHz, A+ = 
2638 MHz, and £i = 200 MHz. 

The entanglement revival behavior occurs even when 
the pair of molecules do not interact, just as a conse- 
quence of the noise due to vacuum electromagnetic fluc- 
tuations. Such fluctuations are always present [28] . and 
the molecules experience a vacuum noise input and radi- 
ate a vacuum noise output. This effect shifts the energy 
levels of the molecules [29 . The so-called Lamb shift has 
been calculated for a specific molecular system in [23 . 
We show such an effect in Fig. [9] where we have consid- 
ered a qubit separation such that V12 = 0, and Ti2 = 0. 
We have introduced a phenomenological A e -shift, start- 
ing from the experimental results reported in [22 , and 
the calculations of Ref. [23] . 

Our result is also in contrast with that obtained in 
Ref. [15 . There, the authors attribute the revival be- 
havior (similar to ours) to the memory (non-Markovian) 
effects induced by two reservoirs acting on the qubits 
in an independent manner. Figure [9] shows that the re- 
vival of entanglement is not, strictly speaking, a non- 
Markovian feature associated to an uncoupled two-qubit 
system. This shows that for the case of independent 
spontaneous emission decay, a similar behavior to that 
reported in [15] can be found. Our result, however, has 
been derived within a Born-Markov description, and as 
such, does not invoke a non-Markovian treatment. On 
the same grounds, this is different from the collapse- 
revival behavior reported for the interacting system of 
Ref. [31 . The behavior showed in Fig. [9] changes with 
increasing emission rates (cf. Fig. [8J: when the sys- 
tem decays more rapidly, the collapse-revival behavior 
disappears. On the other hand, Fig. [8jb) clearly shows 
that, for a fixed rate, collapse-revivals are suppressed for 
higher values of laser intensity, for which ESD becomes 
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the relevant dynamical decay. 




FIG. 10: Comparison between the effects due to the laser and 
the vacuum fluctuations over the entanglement dynamics. In 
the main graph, the solid curve exhibits collapses and revivals; 
£i = 200 MHz, A e = -160 MHz. The long dashed curve 
exhibits ESD; U = 200 MHz, A e = 0. The short dashed 
curve corresponds to £i = 1 MHz, and A e = —160 MHz. The 
initial state is | -00 (1/2)) . The remaining parameters used in 
the simulations are the same as in Fig. [9] By switching to 
the initial state ^/a\00) + y/]5\ll), the inset shows the ESD 
sensitivity to the initial condition; A+ = 20000 MHz, A e = 
— 160 MHz, and £i = 0; the other parameters are as in Fig. [9] 

We have plotted, in Fig. [TUJ the dynamics of the con- 
currence for the initial entangled state \ipo(a = 1/2)), 
to show the interplay between the laser driving and 
the shift experienced by the doubly excited state due 
to background noise fluctuations. The concurrence de- 
cays asymptotically when the laser coupling is absent, in 
agreement with our previous results. 

The main graph of Fig. [To] also shows sudden dis- 
entanglement (but not revival) when the laser coupling 
is strong and there is no contribution arising from the 
vacuum noise (long dashed curve). This behavior is also 
observed in the inset of Fig. [ToJ for a different initial con- 
dition, y^OO) + \fj3\ll): the ESD is more pronounced 
when a tends to zero, and it vanishes for a > 1/2. This 
is because in our density matrix configuration, poo,oo cor- 
responds to the lowest system eigenstate in the standard 
basis, and hence, for a — >• 1, the system has a higher 
probability of decaying to the ground state and the in- 
teresting quantum effects are lost. 

The sensitivity of the ESD dynamics to the initial con- 
ditions is shown in the inset of Fig. [To] In the main fig- 
ure, the entanglement decays asymptotically (uppermost 
curve, £ = 1) for the initial conditions y^Ol) + y/]3\10) 
(for all a), and hence ESD never takes place. However, 



the inset of Fig. 10 shows that for the initial condition 
i/a\00} + y/]3\ll), and absence of laser driving, the en- 
tanglement decays abruptly for values of a close to zero, 
and hence ESD is reinforced at very early times t. This 
clearly indicates that different types of initial entangled 
states experience different input noise from the environ- 
ment. 



An analytical expression for the concurrence, for the 
molecular input values of the inset of Fig. [TUJ follows 
from the calculation of the non trivial density matrix 
elements found for the solution of the master equation: 



Poo.oo = 1 + (1 - a)e~ 2Tt (l - 2e 
Poi.oi = Pio.io = (1 - a)e~ 2Ft (-1 
Pii,ii = (1 -a)e" 2F£ , 
Poo, n = Pn no = Vi 1 ~ a ) a e 



it(27rA e +ir+47rA + ) 



.(6) 



The concurrence C(p) = 2 max {0, |poo,n| — Poi, oil- 
Thus, the ESD time can be directly obtained from the 
condition |poo,n I — Poi,oi < 0; fc> r any a < 1/2, this reads 



1 



1 



a 



a 



(7) 



For a > 1/2 there exists no solution with physical mean- 
ing for t ai in agreement with what is shown in the inset 
of Fig. [To] An experimental study on the engineering 
and modification of the molecules decay rates has already 
been reported [32]; this could facilitate the observation 
of the effects reported in Figs. [9] and [10] 



V. CONCLUSIONS 

We summarize the major findings of the present work 
on dimer entanglement dynamics. We have shown how 
the entanglement evolution, which includes ESD, ESB, 
and collapse-revival dynamics, can be precisely controlled 
by means of an external driving field, and an appropriate 
qubit initialization. In addition, the light-matter inter- 
action can be used to produce transitions between the 
system states in such a way that a reliable production of 
steady state entanglement is favored. There are entan- 
glement dynamics corrections due to the radiative shift of 
the doubly-excited molecular state arising from the vac- 
uum fluctuations: in the scenario where the qubits do not 
experience collective coupling effects (only independent 
interaction with the environment), the entanglement ex- 
hibits a collapse-revival behavior. On the other hand, by 
considering different type of initial conditions, an ESD 
dynamics can also be induced due to the energy shift. 
We have shown that it is crucial, for a correct character- 
ization of entanglement dynamics, to take into account 
the radiative corrections induced by the electromagnetic 
vacuum field, as such shifts or energy corrections perma- 
nently affect the evolution of entanglement. 
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